Incorporating Distributed Debris Thickness in a Glacio-Hydrological Model: Khumbu Himalaya, Nepal

Understanding the future evolution of Himalayan glaciers is important in terms of runoff that provides an essential water source to local populations and has far-reaching downstream impacts. However, the climatic response of glaciers in High-Mountain Asia is complicated by ice stagnation and considerable supraglacial debris coverage, which insulates the ice from warming. Typical runoff modelling only crudely incorporates debris cover and there is currently no consensus on how significantly this may impact future glacier and runoff evolution. Here, a glacio-hydrological model is modified to incorporate fully distributed debris cover, using melt reduction factors that vary depending on debris thickness, and to redistribute mass losses according to observed surface elevation changes. A range of debris thickness data are implemented, including a remote-sensing survey and a modelled debris surface, to analyse the sensitivity of glacier evolution and runoff to possible future debris-cover changes in a series of experiments in the upper Khumbu catchment, Nepal. Simulations are undertaken using climate input data from Regional Climate Model simulations from CORDEX (Coordinated Regional Downscaling Experiment) which are further statistically downscaled using data from the Pyramid meteorological station. Results suggest that the accurate calibration of the model to volume change compensates for the inclusion of distributed debris cover but only if the climatic sensitivity of the calibration period (1999–2010) and the nature of the debris-covered surface remain constant during future simulations. Altering the nature of the debris surface has a significant impact on simulated ice volume, with melt rates under debris suppressed by up to 85 %. The sensitivity of runoff ranges from 60 to 140 million m3 yr-1, although there are considerable uncertainties relating to non-glacial snow melt. Moreover, incorporating locally enhanced melt at ice cliffs into the model also impacts upon volume loss and discharge, with a greater proportion of ice cliffs leading to enhanced volume losses compared to a homogeneous debris surface. Finally, using the most representative model configuration, the future evolution of Khumbu Glacier under various climate scenarios shows continued mass losses with a reduction in volume ranging from 60 % to 97 % by 2100. Runoff trends show an initial increase followed by an eventual decrease, with runoff in 2100 predicted to be 8 % lower than current levels

The hydrology of glacier-bed overdeepenings : sediment transport mechanics, drainage system morphology, and geomorphological implications

Evacuation of basal sediment by subglacial drainage is an important mediator of rates of glacial erosion and glacier flow. Glacial erosion patterns can produce closed basins (i.e., overdeepenings) in glacier beds, thereby introducing adverse bed gradients that are hypothesised to reduce drainage system efficiency and thus favour basal sediment accumulation. To establish how the presence of a terminal overdeepening might mediate seasonal drainage system evolution and glacial sediment export, we measured suspended sediment transport from Findelengletscher, Switzerland during late August and early September 2016. Analyses of these data demonstrate poor hydraulic efficiency of drainage pathways in the terminus region but high sediment availability. Specifically, the rate of increase of sediment concentration with discharge was found to be significantly lower than that anticipated if channelised flow paths were present. Sediment availability to these flow paths was also higher than would be anticipated for discrete bedrock-floored subglacial channels. Our findings indicate that subglacial drainage in the terminal region of Findelengletscher is dominated by distributed flow where entrainment capacity increases only marginally with discharge, but flow has extensive access to an abundant sediment store. This high availability maintains sediment connectivity between the glacial and proglacial realm and means daily sediment yield is unusually high relative to yields exhibited by similar Alpine glaciers. We present a conceptual model illustrating the potential influence of ice-bed morphology on subglacial drainage evolution and sediment evacuation mechanics, patterns and yields, and recommend that bed morphology should be an explicit consideration when monitoring and evaluating glaciated basin sediment export rates

Direct visualisation of collateral ventilation in COPD with hyperpolarised gas MRI

Abstract Background Collateral ventilation has been proposed as a mechanism of compensation of respiratory function in obstructive lung diseases but observations of it in vivo are limited. The assessment of collateral ventilation with an imaging technique might help to gain insight into lung physiology and assist the planning of new bronchoscopic techniques for treating emphysema. Objective To obtain images of delayed ventilation that might be related to collateral ventilation over the period of a single breath-hold in patients with chronic obstructive pulmonary disease (COPD). Methods Time-resolved breath-hold hyperpolarised 3He MRI was used to obtain images of the progressive influx of polarised gas into initially non-ventilated defects. Results A time-series of images showed that 3He moves into lung regions which were initially non-ventilated. Ventilation defects with delayed filling were observed in 8 of the 10 patients scanned. Conclusions A method for direct imaging of delayed ventilation within a single breath-hold has been demonstrated in patients with COPD. Images of what is believed to be collateral ventilation and slow filling of peripheral airspaces due to increased flow resistance are presented. The technique provides 3D whole-lung coverage with sensitivity to regional information, and is non-invasive and non-ionising

Slow dynamics in the 3--D gonihedric model

We study dynamical aspects of three--dimensional gonihedric spins by using Monte--Carlo methods. The interest of this family of models (parametrized by one self-avoidance parameter $\kappa$) lies in their capability to show remarkably slow dynamics and seemingly glassy behaviour below a certain temperature $T_g$ without the need of introducing disorder of any kind. We consider first a hamiltonian that takes into account only a four--spin term ($\kappa=0$), where a first order phase transition is well established. By studying the relaxation properties at low temperatures we confirm that the model exhibits two distinct regimes. For $T_g< T < T_c$, with long lived metastability and a supercooled phase, the approach to equilibrium is well described by a stretched exponential. For $T<T_g$ the dynamics appears to be logarithmic. We provide an accurate determination of $T_g$. We also determine the evolution of particularly long lived configurations. Next, we consider the case $\kappa=1$, where the plaquette term is absent and the gonihedric action consists in a ferromagnetic Ising with fine-tuned next-to-nearest neighbour interactions. This model exhibits a second order phase transition. The consideration of the relaxation time for configurations in the cold phase reveals the presence of slow dynamics and glassy behaviour for any $T< T_c$. Type II aging features are exhibited by this model.Comment: 13 pages, 12 figure

Long-Wavelength Instability in Surface-Tension-Driven Benard Convection

Laboratory studies reveal a deformational instability that leads to a drained region (dry spot) in an initially flat liquid layer (with a free upper surface) heated uniformly from below. This long-wavelength instability supplants hexagonal convection cells as the primary instability in viscous liquid layers that are sufficiently thin or are in microgravity. The instability occurs at a temperature gradient 34% smaller than predicted by linear stability theory. Numerical simulations show a drained region qualitatively similar to that seen in the experiment.Comment: 4 pages. The RevTeX file has a macro allowing various styles. The appropriate style is "mypprint" which is the defaul

Mean flow and spiral defect chaos in Rayleigh-Benard convection

We describe a numerical procedure to construct a modified velocity field that does not have any mean flow. Using this procedure, we present two results. Firstly, we show that, in the absence of mean flow, spiral defect chaos collapses to a stationary pattern comprising textures of stripes with angular bends. The quenched patterns are characterized by mean wavenumbers that approach those uniquely selected by focus-type singularities, which, in the absence of mean flow, lie at the zig-zag instability boundary. The quenched patterns also have larger correlation lengths and are comprised of rolls with less curvature. Secondly, we describe how mean flow can contribute to the commonly observed phenomenon of rolls terminating perpendicularly into lateral walls. We show that, in the absence of mean flow, rolls begin to terminate into lateral walls at an oblique angle. This obliqueness increases with Rayleigh number.Comment: 14 pages, 19 figure

The three-dimensional random field Ising magnet: interfaces, scaling, and the nature of states

The nature of the zero temperature ordering transition in the 3D Gaussian random field Ising magnet is studied numerically, aided by scaling analyses. In the ferromagnetic phase the scaling of the roughness of the domain walls, $w\sim L^\zeta$, is consistent with the theoretical prediction $\zeta = 2/3$. As the randomness is increased through the transition, the probability distribution of the interfacial tension of domain walls scales as for a single second order transition. At the critical point, the fractal dimensions of domain walls and the fractal dimension of the outer surface of spin clusters are investigated: there are at least two distinct physically important fractal dimensions. These dimensions are argued to be related to combinations of the energy scaling exponent, $\theta$, which determines the violation of hyperscaling, the correlation length exponent $\nu$, and the magnetization exponent $\beta$. The value $\beta = 0.017\pm 0.005$ is derived from the magnetization: this estimate is supported by the study of the spin cluster size distribution at criticality. The variation of configurations in the interior of a sample with boundary conditions is consistent with the hypothesis that there is a single transition separating the disordered phase with one ground state from the ordered phase with two ground states. The array of results are shown to be consistent with a scaling picture and a geometric description of the influence of boundary conditions on the spins. The details of the algorithm used and its implementation are also described.Comment: 32 pp., 2 columns, 32 figure

Geomorphological investigation of multiphase glacitectonic composite ridge systems in Svalbard

Some surge-type glaciers on the High-Arctic archipelago of Svalbard have large glacitectonic composite ridge systems at their terrestrial margins. These have formed by rapid glacier advance into proglacial sediments during the active surge phase, creating multicrested moraine complexes. Such complexes can be formed during single surge advances or multiple surges to successively less-extensive positions. The few existing studies of composite ridge systems have largely relied on detailed information on internal structure and sedimentology to reconstruct their formation and links to surge processes. However, natural exposures of internal structure are commonly unavailable, and the creation of artificial exposures is often problematic in fragile Arctic environments. To compensate for these issues, we investigate the potential for reconstructing composite ridge system formation based on geomorphological evidence alone, focusing on clear morphostratigraphic relationships between ridges within the moraine complex and relict meltwater channels/outwash fans. Based on mapping at the margins of Finsterwalderbreen (in Van Keulenfjorden) and Grønfjordbreen (in Grønfjorden), we show that relict meltwater channels that breach outer parts of the composite ridge systems are in most cases truncated upstream within the ridge complex by an inner pushed ridge or ridges at their ice-proximal extents. Our interpretation of this relationship is that the entire composite ridge system is unlikely to have formed during the same glacier advance but is instead the product of multiple advances to successively less-extensive positions, whereby younger ridges are emplaced on the ice-proximal side of older ridges. This indicates that the Finsterwalderbreen composite ridge system has been formed by multiple separate advances, consistent with the cyclicity of surges. Being able to identify the frequency and magnitude of former surges is important as it provides insight into the past behaviour of surge-type glaciers and, if absolute dating is possible, allows for the assessment of surge-type glacier response to climate change on decadal to centennial timescales. Although further investigations into the internal structure of these deposits should be sought where possible, our study demonstrates that geomorphology could be an invaluable tool for reconstructing the formation of composite ridge systems

Distribution and characteristics of overdeepenings beneath the Greenland and Antarctic ice sheets: Implications for overdeepening origin and evolution

Glacier bed overdeepenings are ubiquitous in glacier systems and likely exert significant influence on ice dynamics, subglacial hydrology, and ice stability. Understanding of overdeepening formation and evolution has been hampered by an absence of quantitative empirical studies of their location and morphology, with process insights having been drawn largely from theoretical or numerical studies. To address this shortcoming, we first map the distribution of potential overdeepenings beneath the Antarctic and Greenland ice sheets using a GIS-based algorithm that identifies closed-contours in the bed topography and then describe and analyse the characteristics and metrics of a subset of overdeepenings that pass further quality control criteria. Overdeepenings are found to be widespread, but are particularly associated with areas of topographically laterally constrained ice flow, notably near the ice sheet margins where outlet systems follow deeply incised troughs. Overdeepenings also occur in regions of topographically unconstrained ice flow (for example, beneath the Siple Coast ice streams and on the Greenland continental shelf). Metrics indicate that overdeepening growth is generally allometric and that topographic confinement of ice flow in general enhances overdeepening depth. However, overdeepening depth is skewed towards shallow values – typically 200 to 300 m – indicating that the rate of deepening slows with overdeepening age. This is reflected in a decline in adverse slope steepness with increasing overdeepening planform size. Finally, overdeepening long-profiles are found to support headward quarrying as the primary factor in overdeepening development. These observations support proposed negative feedbacks related to hydrology and sediment transport that stabilise overdeepening growth through sedimentation on the adverse slope but permit continued overdeepening planform enlargement by processes of headward erosion

The Science of Sungrazers, Sunskirters, and Other Near-Sun Comets

This review addresses our current understanding of comets that venture close to the Sun, and are hence exposed to much more extreme conditions than comets that are typically studied from Earth. The extreme solar heating and plasma environments that these objects encounter change many aspects of their behaviour, thus yielding valuable information on both the comets themselves that complements other data we have on primitive solar system bodies, as well as on the near-solar environment which they traverse. We propose clear definitions for these comets: We use the term near-Sun comets to encompass all objects that pass sunward of the perihelion distance of planet Mercury (0.307 AU). Sunskirters are defined as objects that pass within 33 solar radii of the Sun’s centre, equal to half of Mercury’s perihelion distance, and the commonly-used phrase sungrazers to be objects that reach perihelion within 3.45 solar radii, i.e. the fluid Roche limit. Finally, comets with orbits that intersect the solar photosphere are termed sundivers. We summarize past studies of these objects, as well as the instruments and facilities used to study them, including space-based platforms that have led to a recent revolution in the quantity and quality of relevant observations. Relevant comet populations are described, including the Kreutz, Marsden, Kracht, and Meyer groups, near-Sun asteroids, and a brief discussion of their origins. The importance of light curves and the clues they provide on cometary composition are emphasized, together with what information has been gleaned about nucleus parameters, including the sizes and masses of objects and their families, and their tensile strengths. The physical processes occurring at these objects are considered in some detail, including the disruption of nuclei, sublimation, and ionisation, and we consider the mass, momentum, and energy loss of comets in the corona and those that venture to lower altitudes. The different components of comae and tails are described, including dust, neutral and ionised gases, their chemical reactions, and their contributions to the near-Sun environment. Comet-solar wind interactions are discussed, including the use of comets as probes of solar wind and coronal conditions in their vicinities. We address the relevance of work on comets near the Sun to similar objects orbiting other stars, and conclude with a discussion of future directions for the field and the planned ground- and space-based facilities that will allow us to address those science topics